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We summarize some recent ideas relating to anisotropic particle production in high-energy colli¬ 
sions. Anisotropic gluon distributions lead to anisotropies of the single-particle azimuthal distribu¬ 
tion and hence to disconnected contributions to multi-particle cumulants. When these dominate, 
the four-particle elliptic anisotropy C2{4} changes sign. On the other hand, connected diagrams for 
m-particle cumulants are found to quickly saturate with increasing m, a “coherence” quite unlike 
conventional “non-flow” contributions such as decays. Finally, we perform a first exploratory phe¬ 
nomenological analysis in order to estimate the amplitude A of the cos(2(p) anisotropy of the gluon 
distribution at small x, and we provide a qualitative prediction for the elliptic asymmetry from 
three-particle correlations, C2{3}. 
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I. INTRODUCTION 

The observation of large asymmetries, predominantly an elliptic cos(2(/j) asymmetry, in the azimuthal distribution 
of particles produced in heavy-ion collisions has been one of the main indication for the formation of a “nearly perfect 
QCD liquid” [T]. For particles with transverse momenta up to a few times {pr) this phenomenon is usually explained 
in terms of (nearly inviscid) hydrodynamic expansion of an asymmetric “fireball”; on the other hand, for high-pT’ 
particles the asymmetry is thought to originate from energy loss of (mini-) jets along different paths through the hot 
and dense Quark-Gluon plasma [5]. 

More recently, substantial azimuthal asymmetries have also been observed in p-|-Pb collisions at the LHC [SU 
and in d-|-Au collisions at RHIC [3- They are measured via multi-particle angular correlations (see below) and 
were found to extend over a long range in rapidity. By causality, the correlations must originate from the earliest 
times of the collision [5]. The data shows that the asymmetries persist up to rather high transverse momenta, well 
beyond p_l ~ 1 GeV. In fact, a recent publication by the ATLAS collaboration shows that substantial “elliptic” (^ 2 ) 
asymmetries in p-l-Pb collisions at = 5 TeV persist up to p_L = 10 GeV [H]. Final state energy loss is expected to be 
much less prominent in smaller systems created in p-l-p and p-|-A collisions; thus it appears reasonable to investigate 
if azimuthal asymmetries could originate from the instant of collision when (anti-) quarks and gluons are “liberated” 
from the wave functions of the colliding hadrons. Since semi-hard processes involve short-distance QCD dynamics, we 
believe that it is important to develop an understanding of possible origins of azimuthal asymmetries in perturbative 
QCD [IMH]. 

This paper is a write-up of the talks presented by the authors at the “Initial Stages 2014” conference in Napa, 
CA. It is not a comprehensive review but attempts to summarize and combine in one paper a few recent ideas for 
anisotropic particle production and correlations within short-distance, small-a: QCD. 


II. SCATTERING OF A CHARGE OFF A SEMI-CLASSICAL FIELD 


In the eikonal approximation the S-matrix for scattering of a parton in the representation TZ of color-SU(iVc) off 
the target is given by [TS] 

5i(r,b) = ^tr 7 ?,U'''(x) V(y) , (1) 

where r = x — y and b = i(x -|- y) are the dipole radius and the impact parameter respectively. We have implicitly 
assumed that the target field is written in covariant gauge so that the gauge links from x to y and back can be dropped. 
d-ji is the dimension of the representation TZ and V (x) denotes a light-like Wilson line describing the propagation of 
the parton in the field of the target 


U(x) = P exp 



dx A+“(a; ,x)t 



( 2 ) 


Below we shall write most expressions for a fundamental projectile charge, a quark or anti-quark. The S-matrix for 
an adjoint charge (gluon) can be obtained from group theory. 




Nl |5F(r0p - 1 

iV2-l 


( 3 ) 


While S-pif) is complex (for > 3 colors), 5 a( t^ is manifestly real. As we shall see below, this implies that the 
single-particle azimuthal distribution for a quark may in general exhibit odd V 2 n+i moments while that for a gluon 
only has non-zero even moments V 2 n- 

Scattering to high transverse momentum corresponds to small |r|. This allows us to perform a gradient expansion 
of the vector potential A+(a;“,x) resulting in 


5i(r,b)-l = ^^tr (r • E(b)) V ^ 


(i£t.'(r-E(by 


0{r^) 


( 4 ) 


if C-odd exchanges are dropped. The term of order will be used in the computation of C2{3} below but is not 
important for our main point here. In covariant gauge the light-cone electric field of the target in Eq. Q given by 


E\h) 


= / dx-F+^ = - 5 * 


dx A^{x ,b). 


( 5 ) 
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The S-matrix for single parton scattering can be generalized to m particles, 


5m(ri,bi,, 


; 1 — 


jig? 

2Nr 


tr (vi ■ E(bi))^ 


( 6 ) 


where we wrote only the leading order in r to simplify the expression. 

In the current formalism, event averaging corresponds to averaging over the target ensemble, which is defined by 
the field-field correlator. Conventionally, in the McLerran-Venugopalan model [20] one uses 


t(E?(b.)£j(b,)) = 


'(5,,Q2^(bi-b2) 


( 7 ) 


where a general form of the impact parameter dependence of the correlator A(b) with the Fourier image A(k) has 
been introduced. A(b) exhibits a logarithmic divergence as |b| —)■ 0 which is cut off by the dipole scale r since the 
gradient expansion assumes that the electric field is smooth over scales on the order of the size of the probe. 

It should be clear that Eq. Q averages over all fluctuations of the target fields, and hence is isotropic. On the 
other hand, for observables which are sensitive to the angular structure of the target fields, instead we integrate over 
target field ensembles subject to the constraint that the anisotropic contribution to the electric held point in a specihc 
direction a [Hinilll]: 


5 

N, 


|^(£;“(bi)E^(b2))a 


(5“'’Q^A(bi-b2) (%+2^ 


m - 1 


aiQi 



( 8 ) 


That is, we divide the target ensembles into subclasses corresponding to a particular direction of a in the vicinity 
of the point with the coordinates b. The summation over all subclasses (integration with respect to all possible 
orientations a) is performed after the observables (such as the m-particle cumulants) have been computed. In other 
words, the huctuations from one conhguration E(x) to another which spontaneously break 2D rotational symmetry 
constitute slow variables. 

The transverse momentum distribution of scattered partons can now be written as^ 


i2nr 


dN 


kdk d(pk 


= J d-^b J d^re-’^'^-^ S{r,h) 

= J d'^bj dr Sir,ipr,h) . 

The S-matrix satishes 

S{r,(pr) = S*{r,(pr -I-tt) . 

Thus, its real part is even under tpr ^ (fr + (i-C. r —>■ —r) while its imaginary part is odd. 

We can dehne various asymmetry moments of the single-inclusive distribution through 

dN 


Vnihr) = {cosn{(fk - ‘Pa.)) ^ ^ J (^cos{n{pk - Pa)) 


dy kxdkT dpk 


with the normalization 


N = 


dpk / dN 


- 1 / dN 

2tt \ kxdkT dpk / tt \ dk^ 


The brackets (•) indicate an average over all configurations E(x). 

Even (odd) moments have positive (negative) parity under r —)■ —r: 

{cos 2 npk) = + {cos 2 n{(fik + n)) , 
{cos{ 2 n + l)pk) = - {cos{ 2 n + l){(pk + tt)) . 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 


If iS(r, pr) is independent of the orientation of the dipole then all = 0. An angular dependence of its real part gives 
rise to non-zero parity even moments V 2 n', an angular dependence of its imaginary part produces odd moments V 2 n+i- 
Eor a more detailed discussion of the py-dependence of single-particle Vi, V 2 , we refer to Ref. m- We note that 
obtaining non-zero odd-index two-particle cumulants wi{ 2 }, ^ 3 ( 2 } as measured in the experiments is more subtle, see 
sec. IIIIEI below. 


^ Eq. 1^ includes the “no scattering” contribution for transverse momentum exchange k = 0. It plays no role in our subsequent analysis 
since we are interested in finite k only. 
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III. MULTI-PARTICLE CUMULANTS 


A. Connected contributions to high-order cumulants 


We begin this section with the (fully) connected contributions from (Sm) to multi-particle cumulants. We show 
that these generate positive contributions to C 2 {m} and so would lead to complex harmonics V 2 {'m\ if the number 
m of particles is a multiple of four [22JI23]. Hence, that {Sm) gives real V 2 {m} for all m only if the presence of an 
azimuthal anisotropy at the single particle level would generate disconnected contributions. 

Furthermore, we show that |f 2 {m}| beyond m ~ 4 is only weakly dependent on m. This indicates a remarkable 
coherence of the connected “non-flow” contributions obtained from small-a: QCD. Together with their long-range 
correlation in rapidity, the properties are quite unlike “conventional” non-flow, for example, from resonance decays 
or fragmentation of jets. 

We then proceed to discuss contributions from fully disconnected diagrams which arise if rotational symmetry of 
the single-particle distribution is broken. These contribute with opposite sign to C2{2} vs. C2{4}. We also present 
a detailed derivation of the elliptic anisotropy from three-particle correlations, C2{3}. This enables us to analyze a 
“BBGKY-like” hierarchy of m-particle correlations. 

The m-th order cumulant of the elliptic anisotropy is given by 


C 2 {m = 2n} = (exp [i 2{ipi + g)2 - V - <Pn+i - ^^ 71+2 - V2n)\) 


cp. 


(16) 


The normalization in Eq. (161 is dominated by the disconnected contributions, corrections are suppressed by powers 
of 1/iV^. Thus, after averaging with respect to the impact parameters the normalization at leading order in is 


(S'™(ri,...,rm)- I) « ( 


2=1 


(17) 


Equation (|16|) involves all possible contractions that generate the fully connected diagrams. Altogether there are 

-Qr 


(2m — 2)!! contractions: 

(‘S'm (^'l; ; ■ • ■ : bm) 1) 


1 


A(bi - b2)A(b2 - bi) • • • A(b^_i - b„)A(b^ - bi) 

(18) 


4 J (A2 - l)”^-i 
(rir 2 )(r 2 r 3 ) • • • + [(2m - 2)!! - 1] permutations. 


In what follows we adopt a Gaussian A(b) = exp (—b^/cr^) so that 




(19) 


Here 1/Nd is the ratio of the correlated area, 5^, to the area of the projectile, S'_l (the proton in p-A collisions), i.e. 
the inverse number of domains. 

Averaging with respect to the impact parameter and angular variables leads to 


C2{m} = 


m!!(m — 2)!! 


1 


m 2"* [NoiN/ - 1) 


m —1 


The azimuthal harmonics are now readily obtained as [2 

(_I)f-i-i / 1 


(u 2 {m})'" = 


m—1 


with 


mPrn yNoiN^-l) 


/3„ = 2 ^ — 


(m > 2 and even) . 


(m > 2 and even) , 


( 20 ) 


( 21 ) 


( 22 ) 


where jo.fc is the fc-th zero of Bessel function Jo{x). Details on transforming the cumulants, C 2 {m}, to the harmonics, 
U 2 {m}, can be found in Ref. [21]. Eq. (21) also remains true for gluons scattering off the target owing to the cancelation 


of Gasimir factors in normalized observables. 
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FIG. 1: Left panel: absolute value of V 2 {m} for connected graphs only as a function of m. The circles (squares) denote real 
(complex) V 2 {m}. Right panel: V 2 {m} for the disconnected graphs only, see Eq. ( |29[ ). For demonstrational purposes the 
parameter A has been chosen such that the magnitudes of s'-nd V 2 °^{m} are equal at m —>■ oo. Also, we used 

^ = 1/Nn = 1/2. 


The absolute values of the harmonics are approximately equal at large m, quickly approaching the limit 


lim |u 2 {m}| 

m—^oc 


1 jo,l 
Nd 2(iV2 - 1) ■ 


(23) 


Although we were unable to prove this rigorously, we believe that this result holds for any short range correlation 
A(b). The fact that the disconnected “non-flow” small-a; diagrams at large m > 4 are of the “wrong sign” and 
approximately independent of m could perhaps be used to distinguish them from conventional effects. 

The other interesting point here is that the fully connected diagrams give positive cumulants of any order and thus, 
every second V 2 {'m} is complex, starting from m = 4: 


= -C2{4} 


1 

4 


1 



(24) 


This is also illustrated in Fig. A possible resolution consists in an azimuthal anisotropy of the single dipole 
S-matrix [ElEI]. This generates “flow-like” disconnected contributions to the cumulants [22] which we discuss next. 


B. Disconnected contributions to the two- and four-particle cumulants 

Equation ^ corresponds to averaging over all possible configurations of E{b) and is isotropic. However, as we shall 
demonstrate in section EYl for any particular configuration the S-matrix does exhibit an angular dependence, see e.g. 
Fig. In order to account for this anisotropy we instead perform the average according to ([^. 

The first thing to compute is the angular distribution for scattering of a single dipole, for fixed a. Using the leading 
term in Eq. ([^ and Eq. ([^, and performing a Fourier transform to momentum space, as well as an average over the 
impact parameter, one arrives at 


/I dNy^ dN 

\7r dk'^ ) d?k 


l-2A + AA{k- df . 


Consequently, the elliptic harmonic of the single-particle distribution is given by 

V2 = ^ 


(25) 


(26) 
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It is straightforward to generalize the computation of the connected diagrams from above to include the single-particle 
anisotropy. The 2- and 4-particle cumulants turn out to be 


C2{2} 

C2{4} 


Nd ^ 4(1V2 _ 1) 


Nl 


4(iV2 - 1)3 


(27) 

(28) 


The detailed derivation can be found in Ref. [55] . 

Before presenting the result for the 3-particle cumulant C2{3} we first examine the results ( 27|28 ). The first 
term in (27) is the square of the single-particle U 2 ; it is scaled by l/Njj since both particles must scatter from the 
same domain to exhibit a correlation. The second contribution corresponds to genuine non-factorizable two-particle 


correlations, as discussed above. Both contributions are positive; nonetheless Eq. (27) reveals the existence of two 


distinct regimes. For jp the ellipticity is mainly due to the asymmetry of the single-particle distribution induced 

by the E-held domains. In the opposite limit .4 ^ 02 ( 2 } is mainly due to genuine, non-factorizable two-particle 

correlations. 

On the other hand, the fourth order cumulant C2{4} changes sign as a function of A. Furthermore, the magnitude 
of the fully connected contribution relative to U2{1}^ is ~ 1/{A'^N^). Hence, parametrically C2{4} crosses zero when 
A ^ Thus, the presence of both connected and disconnected contributions built from the QCD dipole - 

E-field interaction ~ tr (r • E)^ can in principle describe a change of sign of C2{4} as seen in experiment^. 

We did not manage to derive the general form of C 2 {m} for arbitrary m, if both connected and disconnected 
contributions are included. However, when the single particle contribution dominates, one obtains 


V2{'m} = 




A 

1-1/r; 

D 


(29) 


Consequently, in this case, too, the higher order harmonics are approximately equal to each other 
sufficiently large m. We illustrated this in Fig. a (right). 


V2{m} 


A 

Nd ■ 


for 


C. The three-particle cumulant C2{3} 


In this section we calculate the quadrupole anisotropy from 3-particle correlations 

^^2{3} = C2{3} = (exp + if 2 - ■ 


(30) 


This cumulant is again defined in such a way as to be invariant under a simultaneous rotation of all particle transverse 
momenta by the same angle. 

From Eq. (30) it is clear that the third particle requires a “u 4 -like” structure or else U2{3} would be zero. Such a 
contribution ~ cos(4(^) can be obtained from the expansion of the S-matrix to second order in tr (r • E)^, see Eq. @- 
This leads to the three-dipole S-matrix 


{S 3 ) - 1 = 


1 f (* 5 )' 


(tr (ri • E(bi))2 tr (r 2 • E(b 2 ))^ [tr (rg • E(b 3 ))^]^) . 


2 V ‘2N, 

In this case, the most general decomposition of C2{3} is given by 
C2{3} = (exp {2i{(pi +(f 2 - + (exp (-4i(/?3))(exp {2i{(pi + (^ 2 )) 


(31) 


2 (exp( 2 ji^i))(exp( 2 i(/j 3 ))(exp { 2 i{ip 2 - T’s)))" 
(exp {2i{ipi +ip 2 - 2 (^ 3 )))™""- 


1 -r^P 2 )/ + 2 (exp( 2 f(pi))(exp( 2 i((p 2 - ‘^^Ps)))’' 

(exp {2i{ipi - (^ 3 )))™""' (exp { 2 i{ip 2 - 


(32) 


Although we have computed all of the above terms here we shall focus on the fully disconnected contribution ~ A^ 
as well as on those connected contributions which are of the same order when A = 0{N~^). The second, third and 
the last term in (32) then do not contribute. 


^ Our result probably does not provide a quantitative explanation of the p^-integrated data for C2{4} which is dominated by particles 
with low transverse momenta. Also, the relation of the anisotropy amplitude A and the multiplicity is presently not clear. 
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The overall normalization implicit in (32) will be approximated by the angular average of the fully disconnected 
diagram. It is given by^ 




For the fully disconnected contribution we have 


\ 4 


V 2IVe 




(33) 


C'(a,a')C'(a,a") 


= (1-^ + 2^ (fi • a)^) (1-^ + 2^ (r 2 • a)^) (l - ^ + 2^ (fs • a)^)^ A( 6 i - & 2 )A( 6 i - 63 ) , (34) 

As in ref. [52] here we employed C{a,a') = 27r(5(o — a')A( 6 i — 62 ) and C(a, a") = 2Tr6{a — a")A( 6 i — 63 ) with 
A{bi — bj) = exp(—Averaging over the impact parameters results in 


dH^d^b,dH, ^ ^ ^ 7re7Te_fi 


(35) 


with Nd the number of E-field domains in the target nucleus. 

Multiplying (34) by exp(2z((^i + (p 2 — and averaging over the azimuthal angles leads to the disconnected 

(single particle factorizable) contribution to C 2 { 3 }; in momentum space, 


(exp(2*(<pi + ^2-2<p3))"‘^^- = ^i^^ . 


(36) 


Because we have expanded in the numerator of this cumulant the S-matrix for the third dipole to second order, we 
obtain that C 2 { 3 } ~ drops at high momentum with the square of the px of the third particle. 

The connected and disconnected parts of the fourth term in Eq. (32) are 


Ml 

4iV2 

1 rlrlQi - 


tr (^7^2 ■ E{b 2 )j tr • ^( 63 )) 


A%b 2 - bs) [cOs ((/?2 - P 3 ) + 2A (2 cos (<^2 - Pa) cos {ps - Pa) - cos[p 2 - ¥ 53 ))] 


42 A 2 - 1 

M 

4A2 

^rlrlQt Mh - 63 ) (1 - A - 1 - 2A (fi • a)^) {l- Ap2A {h ■ a)^) 


tr (^i-E{bi)\ \ /tr (fs-Eibs)] \ C{a,a') 

/ d \ / d' 


(37) 


(38) 


respectively. Averaging over impact parameters generates a factor of (1/2A|,). We may now calculate the Fourier 
transform and sum over the 4 contractions of the amplitudes/conjugate amplitudes of the dipoles 1 to 3. This leads 
to 


2(exp(2i(^i))(exp(2iv53))(exp(2i((p2 - T^a)))™™' = ^^^ 2(7/2 _ 1) 


(39) 


The two factors from the fifth term of Eq. (32) each have the form of Eq. (37). Averaging over impact parameters, 
performing the Fourier transform, and summing over the 8 contractions of the amplitudes/conjugate amplitudes of 
the dipoles 1 to 3 leads to 


(exp( 2 i(vi - V^W-- ■ 


Finally, for A ~ 1/Ac we have that 


C2{3} = MM = 


1 IQlfA 


A^ 


Nl4:kM 2 2(A2 - 1) 16(iV2 - 1)2 


(40) 


( 41 ) 


^ In principle one should Fourier transform first to momentum space. At high px this transform is simply F.T.{r} ~ . 





















As already indicated above, we find that C2{3} ^ 1/^1 at high transverse momentum. This is due to the fact that 
in the numerator we expanded the S-matrix to order r| while we only require terms of order in the normalization. 
However, expression (31) for the S-matrix relies again on the gradient expansion of the dipole operator. In sec. IV 


below we shall see that the exact S-matrix (obtained numerically) does appear to include a cos(4(/5) harmonic even at 
order r^, indicating the presence of corrections to the gradient expansion. 

This provides another way for a “?; 4 -like” structure at order r|. In this case C2{3} is given by 


C2{3} = (exp +(p 2 - 2 (^ 3 ))'^*"'=- 

-b (exp +(p 2 - 2 (^ 3 )))““' 


(exp (4z(/93))(exp { 2 i{(pi + 1 ^ 2 ))“™' + 2(exp(2jy)i))(exp ( 2 z(i ^2 - 


(42) 


So far we have not yet computed the diagrams involving a contraction of the third particle with either of the other 
particles. On the other hand, it is easy to write down the contributions from the first two terms in Eq. (42). 


When the third particle is disconnected its S-matrix is given by the S-matrix for a single dipole and we may 
decompose its real part into a Fourier series. 


‘5i(r3) - 1 = A/'(r3) ( 1 -b ^ A2„ cos(2n(</5r - V')) 


(43) 


n—1 


The function Af{r 3 ) = (at small r 3 ) is the isotropic part of the S-matrix and ^ is the “event plane” angle. 

The only term in this series relevant for C2{3} is that for n = 2. The average over (exp(—4 z 9J3)) will contribute 
with 


(exp(-4iv33)) = ^4 y cos(4((/J3 - ipa)) = 


A4 

T 


(44) 


The integrand entering the average over the azimuthal angle for the particles 1 and 2 when they are connected or 
disconnected has the same form as in equations (1^ and (38), respectively. In both diagrams the average over the 


impact parameters will generate a factor of 1/N'^. The overall normalization factor is given by Eq. (33) as before. 
After computing the averages over the azimuthal angles for all particles we have that C2{3} is now given by 


C2{3} 


1 AiA'^ 


(45) 


We assume that A 4 is of order and drop contributions beyond order A'^ ~ If is evident that if the S-matrix 

exhibits a cos(4(^) dependence at order then C2{3} —>■ const at high transverse momentum. Hence, the behavior of 
this cumulant at high px could provide interesting information about the angular structure of the dipole S-matrix. 


D. BBGKY-like hierarchy of m-particle C2 {?ti} cumulants 

In the previous sections we have shown that all C 2 {m} eventually are dominated by the fully disconnected contri¬ 
bution proportional to the single-particle elliptic anisotropy A to the m-th power. This occurs in different stages. 
The four-particle cumulant C2{4} factorizes when (parametrically) A > Nc On the other hand, the two-particle 
cumulant C2{2} requires a stronger E-field anisotropy of order A > N~^. These correlators thus satisfy a BBGKY-like 
hierarchy. On the other hand, in the previous section we have seen that the factorization of C2{3} does not occur at 
some intermediate value of A but, again, for A > N~^, just like for C 2 { 2 }. This correlation function thus represents 
an exception to the hierarchy. In sec. below we shall attempt to go beyond parametric estimates of the connected 
vs. disconnected contributions by performing a phenomenological comparison to data. 


E. Odd-index two-particle cumulants, ci{ 2 } and 03(2} 

In sec. we argued that the angular distribution for a scattered fundamental charge gives rise to odd parity 
moments vi and V 3 . Their dependence on px has been discussed and compared to measured two-particle vi{2}, n3{2} 
in Ref. [H]. We would like to point out here that this issue requires more theoretical investigation, for the following 
reason. 

The two particle correlation function summed over qq, qq, qq and qq channels is C-even if one assumes quark—anti¬ 
quark symmetry of the projectile wave function at small x. Indeed the two-particle S-matrix 

S 2 oc (tr YI(xi) I/(yi) -b trE(xi) Kl(yi)) (trKl(x 2 ) y(y 2 ) -b tr Y(x 2 ) V\y 2 )) 


(46) 
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is real, and so has even cumulants only. 

Therefore, obtaining non-zero ci{2} and C 3 { 2 } may require to account for (at least) one additional soft rescattering 
of the (anti-) quarks besides their hard scattering from the target shockwave. In his talk at this conference Schlichting 
showed that classical Yang-Mills evolution of the liberated gluons in the forward light cone immediately leads to 
non-zero r' 3 { 2 } at time r = 0.1 fm [26j . If such rescattering is soft then the pr-distribution of vi and U 3 shown in 
Ref. |21] should be mostly preserved. Either way, this clearly is an interesting problem which requires more theoretical 
analysis. 


IV. ANISOTROPIC GLUON DISTRIBUTION AT SMALL x 

The main goal of this section is to compute scattering of a dipole off a large nucleus to demonstrate its non¬ 
trivial angular dependence |27j . We shall first consider the classical MV model and then proceed to resum quantum 
fluctuations with large longitudinal phase space via the JIMWLK evolution equation. 


A. Classical McLerran-Venugopalan model 

In the MV model [5^ the large-a: valence partons are viewed as random, recoilless color charges p“(x) described by 
the effective action 




(47) 


with ~ proportional to the thickness of a nucleus; here A denotes the number of nucleons in the nucleus. The 

variance of color charge fluctuations determines the average saturation scale ^ [55]. The Weizsacker-Williams 

fields generated by p“(x) are pure gauges; in covariant gauge. 




(48) 


Using Eq. (48) in Eqs. 0 and 0 we can compute the S-matrix for each configuration of the target fields and extract 
its Fourier harmonics. 

It is rather evident that the random distribution of color charges (x) would generate azimuthally anisotropic soft 
fields. Less trivially, we shall show that the angular structure of the target electric fields does not fluctuate randomly 
on arbitrarily short scales, i.e. that it is characterized by a finite correlation length ~ 1 /Qs in the transverse plane. 
This fact is related to the saturation of the gluon distribution from highly occupied classical fields at momentum 
scales kr < Qs [21]; over distances > 1/Qs the soft, classical color fields become “smooth”. This is why the number 
of domains introduced in previous sections is finite. 

The most crucial aspect, however, is the following. The angular structure will obviously fluctuate from one configu¬ 
ration /9“(x) of valence charges to the next. Averaging over these fluctuations like in Eq. ([^ would obviously project 
onto the isotropic part of the gluon distribution. Instead, we point out that the angular fluctuations of p“(x) are slow 
variables, i.e. that they should be averaged over only after the m-particle cumulants have been computed. 

For a general configuration of the sources, the S-matrix for a fundamental charge is complex. The real (imaginary) 
part corresponds to C-even (C-odd) interactions [55] : 


l-Dp{r) =Re5p(r) = [^^(x) U(y) + U^(y) U(x)] , 

Op(r) = Im5p(r) = ^tr [Ut(x) U(y) - ut(y) y(x)] . 


(49) 

(50) 


We use Monte-Carlo techniques on a lattice describing the longitudinal and transverse coordinates to generate the 
random configurations p“(a:“,x). The number of sites in the longitudinal direction is taken to be = 100, while 
N± = 1024 for either of the transverse directions. We fix the parameters of the lattice such that g'^g.a = 0.05, where 
a = L/N± denotes the transverse lattice spacing. Defining the saturation scale from 


(Sp) (r = V2/Qs) = e- 1/2 (51) 

we determined numerically that Qg « 0.7125g^/r. Further details of the numerical implementation can be found in 
Refs. [miM]. 
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FIG. 2: The configuration-averaged amplitudes (A„}(r) as functions of the dipole size r for n = 1, • ■ • ,4. The fit to (^ 2 ) 
corresponds to the function from Eq. ( |56[ ). Figure from ref. m- 


The azimuthal amplitudes can be extracted by expanding the real and imaginary parts of the S-matrix in a Fourier 
series: 



(52) 


00 


Opir)=JV{r) ^A^„+i(r)cos[(2n-hl)(^^] . 


(53) 


Here, the function Af{r) denotes the isotropic part of the dipole S-matrix. As already mentioned above each amplitude 
contains a random phase which fluctuates from configuration to configuration. To discard this phase we define 
An = f |A(j|; this arises due to 



(54) 


Then, averaging over 10"^ configurations we finally obtain (Ai),-- - , (A 4 ) as well as the variances of Ai and A 2 , 
presented in Fig. 

Our results show that, as expected, the biggest amplitude is the quadrupole; at r < l/Qs the amplitude (A 2 ) ~ 20 %. 
As we argue in the next section, such values are in the range of the asymmetries relevant for phenomenology of high- 
multiplicity p-l-Pb collisions at LHC energies. We stress, however, that in this calculation we did not attempt to bias 
the configurations towards “high multiplicities”, which requires a dedicated investigation. The function (A 2 ) (r) is 
almost independent of r for r < 1/Qs which justifies our treatment in the previous section where A = (A 2 ) has been 
treated as constant. Figure [2 shows furthermore that the variance ^J{{5A2)'^) is similar in magnitude to the mean 
value (A 2 ). This points at rather large fluctuations of A 2 for different configurations. 

Figure shows the same amplitudes as the previous %ure but for E-fields which have been “smeared” over an 
area Trr^ set by the size of the dipole. Comparing Figs, [^and one sees that “smearing” has a negligible effect for 
r < l/Qs while the anisotropy amplitudes at large r are suppressed. This behavior shows the correlation over finite 
transverse distance scales of the angular structure of the E(x) configurations. 

Reference m showed that the MV-model amplitude {A 2 ) (r) matches the distribution of linearly polarized gluons 
(for an unpolarized target) h^^(x,k^) introduced in TMD factorization [32] 



(55) 


Within the framework of the MV model, the result for hf^{x,r) derived analytically in Ref. [32], 



( 56 ) 
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FIG. 3: Same as Fig. for “b-smeared” target E-fields. 


is in good agreement with our numerical results at small values of r 2 Qj^. 

Figure]^ also shows a non-zero amplitude of the cos(4(p) angular component. It appears to be essentially constant 
at small r unlike the ~ behavior expected from the second term in Eq. Q once scaled by M(r) ~ (rQsY at small 
r. This may be due to corrections to the gradient expansion which was used to derive Eq. Q. Such a term would 
provide another contribution to the hexadecupole ti 4 -like asymmetry, c.f. previous section. 

Due to fluctuations of the saturation momentum Qg in impact parameter space |33j every particular configuration 
of semi-classical small-a; fields ( [4^ contains a C-odd component and 0(r) as defined in Eq. (501 is non-zero. This 
results in non-zero odd-index amplitudes Ai and A 3 , see Fig. The figure also shows that the expectation values of 
the odd amplitudes are significantly smaller than A 2 ; as expected, they vanish as r —>■ 0 : 


iO{r) ^ iagT ■ Vb (1 — D{r, b)) ~ I ag Ql Qc cos (pr 


1 - 


/ Ql cos^ pr 


q: 


(57) 


The expression on the r.h.s. corresponds again to a gradient expansion in powers of r, assuming a generic spectrum 
of fluctuations of ( 5 s(b) cut off at [H]. 

The presence of odd harmonics does not indicate that the expectation value of the C-odd part of the S-matrix is 


non-zero. Indeed, the average of the odderon 0(r) over the C-even ensemble generated by the action (47) is zero. 


However, the product of 0(r) with another C-odd operator, which effectively arises due to our dropping of the phases 
of the amplitudes A'^, is even under C-conjugation and its expectation value is not zero. 


B. Quantum fluctuations and high-energy evolution 

In the previous subsection, within the framework of the classical MV model, we showed that scattering of a dipole 
from the soft fields sourced by a particular configuration / 9 “(x) of valence charges is not isotropic, and that the 
amplitudes of the azimuthal anisotropy are quite significant. In this section we consider how these amplitudes are 
affected by small-a; / high energy evolution. This corresponds to a resummation of (nearly) boost-invariant quantum 
fluctuations to the classical field. 

The evolution of the elliptic anisotropy with rapidity was first addressed by Kovner and Lublinsky in Ref. m- 
They solved the BK evolution equation for the dipole scattering amplitude 1 — 5(r) as a function of dipole size and 
orientation. They found that the anisotropy decays exponentially with Y = ln(a;o/a;). Their solution, however, was 
based on the assumption that the impact parameter space is homogeneous. As explained in the previous section, 
even at the level of the initial condition (given by the MV model), the azimuthal anisotropy of 5(r, b) arises due to 
fluctuations of the soft fields in the transverse impact parameter plane. Hence, in this subsection we describe solutions 
of JIMWLK evolution which account for fluctuations of the light-like electric Wilson lines in b-space |27] . 

Going beyond the classical theory, quantum gluon emissions which are enhanced by a large longitudinal phase space 
Y = log Xq/x are resummed by the so-called JIMWLK [331ISS] functional renormalization group evolution. It modifies 
the ensemble of electric Wilson lines over which observables are averaged thereby resumming corrections to all orders 
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in asY. Evolution over a step AY in rapidity opens up phase space for radiation of gluons and modifies the classical 
action (471. The evolution can be formulated in terms of a “random walk” in the space of Wilson lines V (x) 


dYV{x) = V{x) — 

TT 




- - / d^vy(v) 


(x — v)*? 7 *(v) 


(x — u )2 

The Gaussian white noise 77 ® = 77 ® satisfies ( 77 “(x)) = 0 and 

v){y)) = Ois - y). 


(x-v)2 


E1'(v)E(x) . 


(58) 


(59) 


The so-called “left-right symmetric” form of Eq. (58 1 was introduced in Ref. m- We solve Eq. ( [5^ numerically 
assuming a fixed but small coupling as = 0 . 1 ; for such coupling the speed of evolution is at least roughly comparable 
to more realistic running coupling evolution. 

Once an ensemble of Wilson lines on the transverse lattice has been evolved to rapidity E, we can again compute 
the dipole scattering amplitude 5y(r,b), its azimuthal Fourier decomposition and the corresponding saturation scale 
Qs{Y) using Eq. (51). It is important to note here that even though we consider a target of infinite transverse extent 
(periodic boundary conditions), that the evolution equation is solved on a transverse lattice which does allow for 
impact parameter dependent fluctuations. 




FIG. 4: JIMWLK evolution of {A 2 ){r) and {A 4 ,){r) (left) resp. of (Ai)(r) and {A 3 ){r) (right). The lower order harmonics 
correspond to the upper sets of curves. Figure from Ref. [la¬ 
in Fig. 1^ (left) we show the evolution of {A 2 ){r) and (kl 4 )(r) with Y. As already mentioned above, mean-field 
evolution of the dipole was shown to wash out initial elliptic anisotropies rather quickly |12j . On the other hand, here 
we only observe a relatively slow decrease of {A 2 ){r) with Y. This is rather intuitive since both the initial anisotropies 
at E = 0, as well as those of the evolved JIMWLK configurations are generated by fluctuations of the hard “valence 
charges” in the transverse impact parameter plane. Furthermore, we observe that those harmonics which are initially 
small, i.e. (Ai)(r), (A 3 )(r) and (A 4 )(r), in fact increase with E at small r. The harmonics also display universal 
behavior at very large r. The evolution of the amplitudes with E at fixed r Qs{Y) is shown in Fig. 

Thus, we conclude that the anisotropies are not washed out by high energy evolution and that they might be 
essential to describe short distance azimuthal asymmetries observed at LHC. An initial phenomenological analysis is 
presented in the next section. 


V. APPLICATION TO PHENOMENOLOGY OF PROTON-NUCLEUS COLLISIONS 

In this section we present a first phenomenological comparison of the measured U 2 { 2 } and 7 ) 2 { 4 } at high transverse 
momentum to some of the expectations from above. Our analysis is certainly not definitive but preliminary and 
qualitative. Our main goals are: 

• to check if the magnitudes of U 2 { 2 } and 7 ; 2 { 4 } can be reproduced for “reasonable” values of Njy, the number of 
E-field domains, and of A, the E-field cos( 2 ( 73 ) anisotropy amplitude; 
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FIG. 5: JIMWLK evolution of {Ai,2,3a) at fixed rQs{Y). 


• to verify that the connected contributions to the two- and four-particle cumulants indeed describe the splitting 
between V 2 { 2 } and U2{4} observed experimentally at semi-hard px] 

• to estimate the relative magnitudes of connected vs. disconnected contributions to the two-, three-, and four- 
particle cumulants, i.e. how far the respective cumulants are from the factorization limit (dominance of fully 
disconnected diagrams); 

• to make a prediction for U2{3} in pA collisions at the LHC. 


In order to hx A and Nu, we shall use the CMS V 2 {pt) data from 2- and 4-particle correlations in p-|-Pb collision 
at 5 TeV. We focus on the highest multiplicity events. Equations (27) and (28) provide the theoretical expectations 
for the cumulants C 2 { 2 } and C2{4} which we repeat here for convenience: 


c.{2} ^ = r + 

C2{4} ^ -(z;2{4})4 = - 1 


iV3 


(60) 

(61) 


4(iV2-l)3^ 

Note that these expressions do not include subleading corrections in which we defer to a future analysis. Strictly, 




A A 


FIG. 6: Left: Ratio of the full V 2 {m} (m = 2, 3, 4) given by Eqs. ( 27|28|4r I to the disconnected contribution corresponding to 
the first term of each of the equations, respectively. 

Right: Ratio (n2|4|)^^3/(u2|2|)3 as a function of the E-field anisotropy amplitude, A. The dash-dotted line corresponds to 
Eqs. ( [27[^ . The straight horizontal lines are the values from the CMS data (see text for details). 
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Eqs. ( 60|61 ) apply only for A = 0(N~^) and A = 0{Nc respectively. Furthermore, the transverse momenta of 
all particles are assumed to far exceed the saturation scale. 

Depending on the value of A there are two different regimes, see left panel of Fig. for small values of A there 
are strong genuine (non-factorizable) correlations and so the connected diagrams are important; for large values of A, 
however, the cumulants approach the factorization limit where they are dominated by the fully disconnected diagram 
and wher^enuine correlations are suppressed. 

In Fig. (right panel) we plot the A^o-independent ratio (n2{4})^/^/(u2{2})^ as a function of A. The dash-dotted 
line corresponds to Fqs. (|^ 28) while the straight horizontal lines represent that same ratio obtained from the two 
highest pt data points for vAAi '*^2{4} shown in Fig. As one can see the high-pT data allows two regimes of 
A: one around A ~ 0.2 and another for 0.35 ^ ^ ^ 0.7. From Fig. ^we see that the first solution is in the regime 
where strong correlation effects are present while the second one is close to the factorization limit. 

The comparison of Fqs. ( 27|28 l to the CMS data in the high px region is shown in Fig. The values of A — 0.20 
and A = 0.53 employed in the figure correspond to the two possible solutions mentioned in the previous paragraph. 
Njj was fixed so as to reproduce the correct magnitudes of ^ 2 ( 2 } and U2{4}. Since both set of parameters, for small 
and large A, are able to describe the data with comparable quality we must conclude that our analysis is not sufficient 
to determine A and Nd uniquely from the high px data alone. If the data down to about px = 1 GeV is included in 
the analysis then the model would prefer smaller values A ~ 0.2 [2T| . 
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FIG. 7: V 2 {pt) from two- and four-particle cumulants for A = 0.20 (left) and A — 0.53 (right), respectively. In either case Nn 
was fixed in order to reproduce the magnitudes of U2{2} and U2{4} simultaneously. The splitting between these two quantities 
is due to the contribution from connected diagrams. Data by CMS [H]. 

Figureshows the predictions for U2{3} obtained from both expressions 


1 I Ql fA 


c,{3} = (.,{3})^^_-:^ V 


iV2 4fc| 


A^ 


2 2{N^ - I) 16(Af2 - 1)2 


C2{3} = 


1 AiA"^ 


(62) 

(63) 


derived above, using the same values for A and Nd as deduced from U 2 I 2 } and t'2{4}. Recall that ( [6^ results from 
an expansion of the dipole S-matrix to second order in tr (r • E)^ while (63) arises if the S-matrix exhibits a cos(4(^) 
component already at order A. For Fig. |^we assumed that A 4 = A^ to avoid introducing additional parameters. The 
figure shows that despite the uncertainty in A and No that 112 {3} does not vary too widely. We expect U2{3} « 2 — 4% 
for semi-hard transverse momenta. 


VI. SUMMARY AND OUTLOOK 

Our goal here was to provide a summary and overview of recent ideas regarding anisotropic particle production at 
semi-hard transverse momenta in high-energy collisions. The basic point is that azimuthally anisotropic correlations 
should occur essentially due to an anisotropic small-T gluon distribution. 

The McLerran-Venugopalan model for the gluon distribution of dense hadrons or nuclei integrates out the fast 
dynamics of the large-x degrees of freedom and replaces them by “frozen” sources for the small-a; semi-classical fields. 
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We point out that each such configuration exhibits azimuthal anisotropies with a finite transverse correlation length, 
and that the angular structure of these configurations is a slow variable, too. 

As a consequence, before one averages over the random angular structure of the source, the single-particle dis¬ 
tribution due to scattering of a projectile parton off such a target is anisotropic. This gives rise to contributions 
to multi-particle correlations from disconnected diagrams [T2 HU |22]. By analogy to the BBGKY hierarchy the 
disconnected contributions dominate the m-particle correlation functions in the limit of large anisotropy A of the 
gluon distribution of the target. More specihcally, they have been shown [52] to lead to a sign flip of the four-particle 
“elliptic” cumulant C2{4}. 

On the other hand, the connected contributions to the cumulants from small-cc dynamics exhibit a rather unex¬ 
pected^ coherence in that C 2 {m} depends weakly on the order m of the cumulant, for sufficiently large m |23j . quickly 
approaching a constant as 1/m —>■ 0. Also, unlike conventional “non-flow” contributions, the connected diagrams 
from the CGC (coherent small-x QCD dynamics) are long range in rapidity [51 ITT]. 

Much work is still needed before we might claim to understand the data. From the point of view of phenomenology 
one should, for example, compute subleading in N~'^ and Qg/Px corrections to the cumulants in order to improve 
the analysis performed in sec. [V] It would be very interesting, too, to resum the time evolution of the gluon fields in 
the forward light cone [26] in order to understand over what range of px final-state interactions are important. Also, 
this approach could clarify the time scale over which odd-index cumulants like ci{m} and C 3 {m} develop. As a last 
point, let us mention that the effect of the multiplicity bias employed in experiments on the anisotropy of the gluon 
field configurations is poorly understood at present. 

An interesting theoretical issue is the relationship of the azimuthal cumulants to the gluon distributions introduced 
(via specific operator-level relations) in TMD factorization [31]. Within the MV model at least, it has been shown 
by explicit computation that M(r), which corresponds to U 2 { 1 }(pt) A high px, coincides with the distribution of 
linearly polarized gluons hf^{r) [55]. It remains to be seen if this relation still applies after resummation of small-x 
quantum corrections. Azimuthal correlations in high-energy p-l-p and p-|-A collisions could provide much insight into 
non-trivial QGD dynamics. 
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